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We study the transverse momentum spectra of fermions and bosons produced in strong, time- 
dependent abelian field. The transverse size of the abelian field is finite, similarly to color strings and 
ropes. Difi'erent time-dependent field strengths are investigated in a kinetic model, and transverse 
momentum spectra are calculated for fermions and bosons. These spectra display exponential or 
polynomial behavior at high pr, depending on the given time dependence. We compare our spectra 
to lattice result for a classical gluon field and obtain surprisingly good agreement in certain cases. 
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I. INTRODUCTION 

Ultrarelativistic heavy ion collisions at SPS and RHIC 
energies (i/s = 10 — 200 AGeV) provided an enormous 
volume of experimental data on particle production 0. 
The microscopic mechanisms of hadron production are 
not fully understood and many different descriptions co- 
exist to explain these data. Transverse momentum spec- 
tra of produced hadrons have been measured in a wide 
momentum region (e.g. < pT < 15 GeV at RHIC), and 
can become a decisive test between different perturbative 
and non-perturbative models of hadron formation. 

Investigations of pp collisions at high energies led to the 
introduction of chromoelectric flux tube ('string') mod- 
els, where these tubes are connecting the quark and di- 
quark constituents of colliding protons Q, y, Q . Strings 
are good examples of how to convert the kinetic en- 
ergy of a collision into field energy. When these flux 
tubes become unstable, new hadrons will be produced 
via quark-antiquark and diquark-antidiquark pair pro- 
duction. Such models can describe experimental data 
successfully at low px, Pt < 2 — 3 GeV. At higher px 
perturbative QCD-based models are working 0,|a,l3|- 

In heavy ion reactions finite number of nucleons collide 
and the number of produced strings scales with the num- 
ber of participant nucleons. Since hadron production at 
low pt scales with participant nucleons in a wide energy 
range, string models could reproduce available data sur- 
prisingly well in the soft region at SPS energies 
However, the obtained high density for strings strength- 
ened the idea of string melting and the formation of color 
ropes 1^, which fill the finite transverse extension par- 
tially or completely. Following these ideas, measured 
strangeness enhancement was explained successfully by 
rope formation 4] . This result has indicated the possible 
formation of extended color fields at SPS energies. 

At RHIC and LHC energies the string density is ex- 



pected to be so large that a strong gluon field will be 
formed in the whole available transverse volume. Fur- 
thermore, the gluon number will be so high that a clas- 
sical gluon field as the expectation value of the quan- 
tum field can be considered and investigated in the reac- 
tion volume. The properties of such non-abelian classical 
fields and details of gluon production were studied very 
intensivelyduring the last years, especially asymptotic 
solutio ns M [Toll . Fermion production was calculated re- 
cently |llLll2j |. Lattice calculations were performed also 
to describe strong classical fields under finite space-time 
conditions p^lT^. 

Fermion pair production together with boson pair pro- 
duction were investigated by kinetic models of particle 
production from strong abelian J^5, 16, 17, 18, 19, 2^, 
I21L l23 | and non-abelian fields. These calculations 
concentrated mostly on bulk properties, the time depen- 
dence of energy and particle number densities. 

Our main interest is the transverse momentum distri- 
bution of produced fermions and bosons. Before perform- 
ing non-abelian kinetic model calculation, we would like 
to understand the role of time dependence, the interplay 
between production and annihilation rates in a kinetic 
approach and the influence of finite transverse size on 
the transverse momentum distributions. 

In this paper we summarize our results applying a 
kinetic model with a time-dependent abelian external 
field characterized by finite transverse geometry. We 
concentrate on transverse momentum spectra for pro- 
duced particles. Section 2 summarizes the field theo- 
retical backgrounds for boson and fermion production in 
strong abelian field. The kinetic equation is explained 
briefiy in Section 3. In Section 4 the influence of time 
dependence on fermionic and bosonic transverse momen- 
tum spectra is displayed and the obtained results are 
compared to lattice calculations. In Section 5 we discuss 
our results. 



2 



II. BOSONS AND FERMIONS IN EXTERNAL 
FIELD 

Let us consider a massive boson field (j) in an external 
classical abelian vector field, A^^. The Lagrangian 



The energy of the bosonic quasiparticles reads 



C^ = Dlcp*D''(l)-m\cp*cp 
leads to the equation of motion 



(1) 



(2) 



where Z?^ = + ie^A^ with bosonic charge e^. The 
bosonic mass is labelled by m-(-. 

We will choose a longitudinally dominant vector field 
in Hamilton gauge with the 4-potential = (0, 0, 0, A^), 
which is the most appropriate for our investigation in the 
ultrarelativistic energy region. 

To imitate confinement properties of Yang-Mills fields, 
the component A^ is limited in the transverse direction, 
and a finite 'flux tube' radius tq is considered. The ex- 
ternal field is cylindrically symmetric. It vanishes out- 
side the tube, A^{t,r > rQ,ip,xz) — 0, and it is ho- 
mogeneous inside the flux tube, A^{t,r < VQ^ip^x^) = 

AHt) Eiiiiiii. 

The Klein-Gordon equation Q reads for the boson 
field r, (/?, X3) as 



- Vi - di - 2ie+A3 {t)d3+elAl (t) +m^_^ 



= 0, (3) 



where the transverse Laplace operator is given by 
92 Id I 



(4) 



For the bosonic field function in eq. we are looking for 
the following solution: 



T(t)R{i 



(5) 



One equation is obtained for the time-dependent part. 



+ (fc3-e+A3(t))Ve' T(i) = 0, (6) 



and one for the spatial dependence, 

r^R\r) + rR{r) + (eV^ - t^) R{r) 



0. 



(7) 



Here e is the separation constant which will be fixed later. 

Considering flux-tube boundary condition for A^it) 
and the constraint R{erQ) = on the surface of the flux 
tube, the solution for the boson field is 



^/nroJi (e„;ro) 



where e„;ro is the n*'' zero of the Bessel function J/ and 
the constant J; (e„;ro) appears during the normalization 
of the field function. 



uj1{t,kT,k3) = 



p2 



(9) 



where P+ = k^, — e+A'i{t) is the kinetic longitudinal mo- 
menta and e„i labels the discrete transverse momenta. Iw 
Section 3 the transverse momenta spectra of the newly 
produced bosons will be determined at these discrete 
transverse momentum values, e„j. 

Massive fermions (77i_,e_) can be described similarly 
in the presence of the above external classical abelian 
field Considering the fcrmionic Lagrangian 



one obtains the corresponding equation of motion 

{i-f^'^f, - e^-f^Asit) - m^)iP{x) = . 
The solution of eq. is wanted in the form of 



(10) 



(11) 



(12) 



where an auxiliary field ip{x) was introduced |2J|. In the 
cylindrically symmetric case, one obtains the following 
equation for ip{t, r, 0, 2:3): 

[di - Vi - dl - 2^e_A,{t)^3 + elAl{t) 

+ie--/°-f^doA3 + = . (13) 

Now we can look for the solution of ea. H13|) in a form 
similar to eq. : 

^{t, r, if, xs) = Ne'''^e-''''^'f{t)R{r)xx ■ (14) 

Eigenvector of matrix 7073 is labelled by xa, A = ±1. 
During the solution of ea. (|13|l . the radial equation for 
R{r) leads to eq.©. Thus, one obtains R{r) = R{r) 
and the same radial solutions with the same separation 
constants e„; as above. The one-particle energy for the 
auxiliary field is the following: 



QUt,kT,k3) = ml + ell + Pi 



(15) 



Here P-{t) = k^ — 6-^3 (t) labels longitudinal momenta. 
Time-dependent component T{t) satisfies the equation 

f{t) + (t, fcr, fc3) + ie^Aait)] f{t) = . (16) 

Using ea. H12l) for the solution of the auxiliary field, 
one can reconstruct the fcrmion field Tp(t,r,ip,X3). The 
dispersion relation in ca. (|15|l remains valid for these 
fermions. 

The finite radius of the flux tube determines a minimal 
transverse energy, eio — 0.2404/ro. For pp collisions we 
obtain eio(ro = 1 fm) = 47.4 MeV^, for central AuAu 
collisions it is eio{ro = 6.5 fm) = 7.3 MeV^, which are 
very small, but nonzero values. 
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III. THE KINETIC EQUATION 

The kinetic equation with time-dependent external 
field can be constructed on the basis of the calculated 
bosonic and fermionic field functions. Using a time- 
dependent Bogoliubov transformation approach (see e.g. 
Refs. I24. l25| ') or a quasiparticle representation ap- 
proach |2ll |. one obtains the following kinetic equation: 



df±{t,k3 
dt 



e±E{t) 



d.f±{t,h) 
dks 



S±{t,h). (17) 



Here E{t) — ~A3{t). One-particle distribution functions 
for bosons(-l-) and fermions(— ) are labelled by f±{t,k^). 
We note that the measurable longitudinal momentum of 
the produced particle has the value of P± . 

Equation (|17|) will be solved at the discrete transverse 
momenta e„i, which depend on the transverse radius (rg) 
of the chromoelectric field. Indeed, eq. (|17|l can be solved 
at continuous transverse momenta and the solutions at 
corresponding e„; will be selected. 

The source term S±{t^ k^) contains information about 
the time evolution of the system and is determined as 



S±{t,k^)='^W±{t,k:i) I dt' W±it,t';k3)- 



[l±2/±(i',fc3)]cos 



where 
ujl{t,t';k3) 

P±{t,t') 
W±it,t';k3) 



2 J dr uj±{t, r; k^) 
—1 



= m| + e2^+P±(t,t') 
t 

= + e± y dTE{T) , 
t' 

^ e±E(t)P±{t, t' 
^±{t,t'-ki) 



P±{t,t') 



(18) 

(19) 
(20) 

b 

(21) 



Overlined notation indicates appropriate asymptotic 
behavior, namely W±{t, k^) = W±{t, t' —00; fca). Sim- 
ilar asymptotics are valid for P±{t) and u}±{t). 

The difference in boson and fermion production ap- 
pears through the spin factor in eq. I|21|l . namely 5+ = 
for bosons and s_ = 1/2 for fermions. 

Now we introduce new auxiliary functions, v± and u±, 
and rewrite the integro-differential eq. H17() in a more con- 
venient form containing full derivatives: 



f± 
v± 
u± 



^W±v± , 

T^±(l±2/±) 
2uj±v± . 



2u!±u± 



(22) 

(23) 
(24) 



This system of differential equations is solved for all sets 
of {e„i, fcs}, which correspond to the different transverse 
and longitudinal momenta. 



The functions f±,v±,u± start with zero initial value, 
describing a system without any bosons or fermions at 
the beginning. Assuming different time dependences for 
the external field, E(t), we obtain different particle pro- 
ductions. We would like to note that produced charged 
particles generate an inner field which interferes with the 
external one This effect is neglected here, as well 

as a possible collision term, C±(f, fcs), on the right-hand 
side of eq. (|17(l '2^ . A given external field can mimic a 
certain self-consistent solution as expansion, collision or 
back reaction infiuenced time dependence of the field. 



IV. NUMERICAL RESULTS 

In heavy ion collisions, one can assume three differ- 
ent types of time dependence for the chromoelectric field 
to be formed: a) pulse-like field develops with a fast in- 
crease, which is followed by a fast fall in the field strength; 
b) formation of a constant field (Eq) is maintained af- 
ter the fast increase in the initial time period; c) scaled 
decrease of the field strength appears, which is caused 
by particle production and/or transverse expansion, and 
the decrease is elongated in time much further than the 
pulse-like assumption. 

Figure ^displays three sets for the time dependence of 
the external field to investigate the above three different 
physical scenarios numerically: 



Epuiseit) ^Eo ■ [1 - tanh(t/<5)] 

Epnlse (t) at t <0 

Eq at t>0 



Econst {t) 
Escaled{t) 



Epulse{t) Sit t <0 

JT+tftoF a* * ^ 



(25) 
(26) 

(27) 




FIG. 1: The time dependence of external field E{t) in three 
physical scenarios: a) pulse {dotted line); b) constant field, Eq 
{dastied line); c) scaled decrease {solid line). 
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FIG. 2; The transverse momentum spectra for bosons in the 
three physical scenarios (see Fig. and text for explanation). 
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FIG. 4: Longitudinal momentum spectra for fermions (thick 
line) and bosons (thin line) at kr/Ey^ — 1 and t — 



1/2 

Here we choose 5 = 0.1/ Ef^' , which corresponds to 
RHIC energies, and k = 2/3 to indicate a longitudinally 

1/2 

scaled Bjorken expansion with to = 0.01/ Eq' . In fact, 
the whole time dependence is scaled hy Eq ' . 

Figures |21 and O display the obtained transverse mo- 

1 /2 

mentum spectra in midrapidity at time t = 2/Eq for 
bosons and fermions. Pulse-type time dependence leads 
to exponential spectra {dotted lines), f± cx exp(fcT/7±), 

1/2 

where the slope value is T+ = 1.6 • ii^g for bosons and 

1 /2 

T_ = 1.45 • Eq' for fermions. In the other two cases, we 
obtain non-exponential spectra, which are generated by 
the long-lived electric field. The transverse spectra from 
constant {dashed lines) and scaled {solid lines) fields are 




kr/E,' 



FIG. 3: The transverse momentum spectra for fermions in the 
three physical scenarios (see Fig. and text for explanation). 



close, because the production and annihilation rates bal- 
ance each other. Slight differences appear because of the 
fast fall of the scaled field immediately after t = 0. We 
note here, that indeed we have discrete transverse mo- 
menta, e„i, but with such small steps, which allow us to 
display continuous transverse momentum spectra. 

Comparing the boson and fermion production in 
midrapidity, one can recognize the large difference in the 
production rates and see a clear fermion dominance at 
all transverse momenta. The reason of this dominance 
can be understood from Fig. 4, which displays the longi- 
tudinal particle spectra, f±{P±), at fixed transverse mo- 
menta, kx/Ef/^ = 1, and time, t = 2/Eq^'^. Fermion 
spectra (thick line) have a reasonable value around zero 
longitudinal momenta (which means midrapidity), but 
boson spectra (thin line) display a 'dip'-like structure 
here. The presence of this 'dip' is related to the struc- 
ture of the transition rate W{t, t'; k^) in eq. (|?T|l and the 
longitudinal dominance of the electric field defined via 
A3{t): at zero longitudinal momenta {P±) there is no 
boson production, but fermion production is character- 
ized by a well defined finite value. 

At extremely high-energy heavy ion collisions we ex- 
pect this longitudinal dominance of the chromoelectric 
field. If transverse electric field components become com- 
parable to the longitudinal one at lower energies, then the 
boson yield will be much larger. The presence of scatter- 
ing term, C±, can increase the yield in midrapidity from 
the neighbouring rapidity cells, where the yield is large, 
even if a longitudinal field component, , exists alone. 

In Refs. |il3j transverse momentum spectra for giuons 
were determined in SU(3) classical field theory. Com- 
parison between this non-abelian and our abelian calcu- 
lation can be made, if we fix both field strength at the 
same value. The energy density of the classical gluon 
field is given in Refs. |l3 as esu3 = O.lTA^/g^, where 
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FIG. 5: Transverse momentum spectra for fermions {thick 
lines) and bosons {thin lines) from our calculation with scaled 
time evolution, Escaied{'t), and gluon production {diamonds) 
from SU(3) lattice-QCD calculation [l^. 

g — 2 {as — 0.33). In our calculation, the initial en- 
ergy density is related to the field strength Eo{t = 0): 
e — E'q/2. Assuming the same 'initial' energy den- 
sity in the abelian and the non-abelian case, we obtain 

1 /2 

E^' /As = 0.54. Now we can rescale our transverse mo- 
mentum spectra from Figures |2] and 13 - we choose the 
results with scaled time-dependent field, Escaied (solid 
lines). In Figure |31 our results at RHIC energies for 
fermion {thick solid line) and boson [thin solid line) are 
displayed together with non-abelian result for classical 
gluons {open diamonds). Our fermion spectra and 
the lattice result on classical gluons overlap in the mo- 



mentum region 2 < kr/As < 10. Possible parametriza- 
tion of the two spectra are very close to each other and 
display the \og{kT / As) ■ [Ut / As)~'^ behavior known from 
Ref. 1^ and from perturbative QCD "271. An accept- 
able slight modification of the time evolution in eq. H27() 
could lead to a full overlap, even at high-fcr. This can 
be demonstrated at kr/As > 15 with S ~ 0.01/ Eq for 
fermion production {thick dashed line). 



V. DISCUSSION 

In this paper we investigated fermion and boson pro- 
duction from a strong abelian field in a kinetic model. As- 
suming ultrarelativistic limit and a longitudinally domi- 
nant external field, boson production remains relatively 
weak in the mid-rapidity region. However, fermion pro- 
duction is not effected by this setup and a fermion domi- 
nance appears in our kinetic model. Introducing different 
parametrizations for the time evolution of the external 
field, we obtain different transverse momentum spectra 
for fermions and bosons. Only time evolution determines 
if these spectra are exponential or polynomial. Investi- 
gating fermion transverse momentum spectra, our results 
with abelian external field overlap with the spectra ob- 
tained in lattice-QCD calculations for classical gluons at 
the same energy density. 
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